arXiv:1504.06306v2 [hep-th] 25Jun2015 


Preprint typeset in JHEP style - HYPER VERSION 


On Multi-step BCFW Recursion Relations 


Bo Feng“^, Junjie Rao“, Kang Zhou“* 

Zhejiang Institute of Modern Physies, Zhejiang University, Hangzhou, 310027, P. R. China 
’’Center of Mathematical Science, Zhejiang University, Hangzhou, 310027, P. R. China 


Abstract: In this paper, we extensively investigate the new algorithm known as the multi-step BCFW 
recursion relations. Many interesting mathematical properties are found and understanding these aspects, 
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the correct answer, without recourse to any specific knowledge from field theories, besides mass dimension 
and helicities. This process consists of the pole concentration and inconsistency elimination. Terms that 
survive inconsistency elimination cannot be determined by the new algorithm. They include polynomials 
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the Standard Model plus gravity to illustrate its power and limitation. Ensuring its workability, we also 
tentatively discuss how to improve its efficiency by reducing the steps. 
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1. Introduction 


In the past decade, the BCFW recursion relation [1, 2] had been an efficient on-shell method to calculate 
tree-level scattering amplitudes. Pedagogical reviews on this topic can be found in [3, 4], Still, it encounters 
certain difficulties when there exists no ‘good’ deformation as those found in [5, 6], i.e., the real amplitude 
does not vanish under the large z limit, where z is the deformation parameter. The recursion relation then 
fails to capture a residual part called the boundary term, which corresponds to the residue at infinity of 
the deformed amplitude. 

Many related studies have been achieved including: introducing auxiliary fields to eliminate boundary 
terms [7, 8], analyzing Feynman diagrams to isolate boundary terms [9, 10, 11], expressing boundary terms 
as roots of amplitudes [12, 13, 14], collecting factorization limits to interpolate boundary terms [15] and 
using other deformations for better large z behavior [16]. 

Recently, a new algorithm named as the multi-step BCFW recursion relations [17] was established to 
tackle this problem universally. Its major idea of using auxiliary deformations can be traced back to [18], 
while the latter aims for one-loop amplitudes. This approach considerably widens the category of quantum 
field theories of solvable tree amplitudes by using BCFW deformations only [19]. However, some common 
puzzles encountered in practice still lacks a formal study. One core question is: How to reach the correct 
answer within finite, definite steps, if an amplitude is solvable by the algorithm? 

In this paper, we will first explore multi-step BCFW recursion relations by investigating the algebra 
of BCFW deformation generators and the commutativity of constant extractions. Next, we will seek for a 
universal approach to reach the answer and ensure that it is correct. This safety promise relies on very little 
knowledge of a particular QFT, besides mass dimension and helicities, hence the algorithm is expected to 
be able to solve for all massless tree amplitudes, with certain limitation as addressed below. 

It is well known that on-shell methods heavily rely on factorization properties of amplitudes, and the 
latter is a reflection of locality and unitarity. These properties are mathematically implemented on poles 
of amplitudes and their residues. For amplitudes that admit polynomials, no on-shell methods so far can 
fix this ambiguity. One can list all possible forms of polynomials as basis, but to determine the coefficients 
will unfavorably call for more traditional means such as Feynman rules. In this work, we will clarify the 
applicable range of multi-step BCFW recursion relations and explore all possible forms of polynomials and 
their generalized cousins called pseudo polynomials and saturated fractions. The latter two objects can be 
fixed by other types of deformations, and having them fully identified is in fact useful. 

The paper is organized as follows. In section 2, we review the multi-step BCFW recursion relations 
and explore the commutativity of constant extractions. In section 3, we propose the systematic process to 
calculate amplitudes after finite, definite steps, and clarify its applicable range and limitation. In section 
4, we apply it to the (massless) Standard Model plus gravity to demonstrate its workability. 

2. Multi-step BCFW Recursion Relations 

In this section, we briefly review the multi-step BCFW recursion relations, in the novel language of extrac- 
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tion operators. After that, the commutativity of constant extractions will be explored. 


2.1 Extraction operators 

For a general BCFW deformation {ai\bi] (only two legs are shifted for each i), namely 


Aai ^ai — ZiXbi, Aft. —)• Aft. + Zi\ 


ai ) 


( 2 . 1 ) 


let’s define two operations on an amplitude-like rational function i?(Ai,Aj) via^ 

^^4^] = - E 


finite 


dz' r dz' " ~ 

- R{Xai — ZiXb^, Xbi ZiXai), Ci[R] = ® - R{Xai — ZiXb^, Xb^ + ZiXa^), (2-2) 


Zi 


Zi 


where Pi and Ci are the pole and constant ‘extraction operators’, which capture residues at finite locations 
except zero and infinity respectively. For a real amplitude A, Pi can capture its physical poles only. But 
a general R, such as Pi A = Pi[A] or CiA = ^[A], may also contain spurious poles, which is well known. 
Therefore the detectable poles (those which have dependence on Zi, as defined in [17]) at finite locations 
can be either physical or spurious. 

By definition Pi + Ci = /, where I is the identity operator. When we calculate an amplitude, starting 
by the 0th step, the amplitude is unknown, so is the Cq operation. However, the Pq operation represents 
exactly the BCFW recursion relation, hence we actually reconstruct this part by employing factorization 
properties, rather than manipulating the unknown amplitude. Conventionally, Co is called the boundary 
term with respect to Pq, which will be dissected into many parts to be determined. The dissection means, 
by expanding I for (n -|- 1) times repeatedly, we have 


I — Pn CnPn-l + . . . + CnCn-1 ' ' ' C 2 P 1 + CnCn-1 ‘ ‘ ‘ C'2C'iPo + C'nC'n-1 ‘ ‘ ‘ C 2 C 1 CQ, (2-3) 


note that I always acts on A implicitly. If the final boundary term CnCn-i ■ ■ ■ C 2 C 1 C 0 vanishes, we have 


I — Pn p CnPn-l -|- . . . + CnCn—l ' ' ' C 2 P 1 + CnCn—1 ' ' ' C 2 C 1 P 0 . (2-4) 

This identity formally represents the ‘multi-step BCFW recursion relations’. Importantly, the workability 
of this multi-step approach relies on the existence of a sequence of deformations numbered by 0, 1 ,... ,n 
for which CnCn-i ■ ■ ■ C 2 C 1 C 0 = 0. The latter is the key condition we will mainly focus on. 

The operators above have a general algebraic property, namely the projectivity: 

CiC^ = a, (2.5) 

to prove this, we first explicitly expand the deformed R as^ 

^(-•) = E 7 -+ ( 2 . 6 ) 

k (®0fc + 0,lkZi + Cl2kZi ) I 

^Pi and Ci used here are identical to V- and C- in the appendix of [17]. 

^In practice, one can use the ‘Apart’ function in Mathematica to separate the pole and regular terms with respect to z. 
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with dfc > 1. In the expansion, when a 2 k vanishes, hik must also vanish, otherwise a linear recombination 
of the numerator can further lower by one^. Now observe that performing the same deformation twice 
is equivalent to replacing Zi by (zj + z'), as 


.') = R(z, + = V-^ 

k (“Ofc + + z[) + a2k{Zi + z[YY’‘ 


cq+ ''^^^ci{zi + z[)\ (2-7) 


then 


^ / —R{zi,z[)= / ^ [co + ^czz'M =co. 


OO '^i J OO '^1 


hence CiCiR = CiR = cq. By using Pi = I — Ci \i is trivial to find that 

PiPi = Pi, CiPi = PiCi = 0 . 

Besides projectivity, a more intricate property is the commutativity: 


( 2 . 8 ) 


(2.9) 


CiCj = CjCi, (2.10) 

which demands certain condition, as will be investigated shortly. If it holds, again with Pi = I — Ci one 
can hnd that 

PiPj = PjPh CiPj = PjCi- ( 2 . 11 ) 

When all C’s are chosen to commute with each other in the expansion (2.4) for a particular amplitude, 
each term is ‘orthogonal’ to the others. This orthogonality has a nice meaning: Each term contains non¬ 
overlapping pole terms, consequently one can capture all pole terms step by step without checking whether 
the previous parts are disturbed by new operations. While commutativity may considerably simplify the 
calculation, it is obviously not necessary for (2.4) to work. 

One last digression is when we do practical calculations, it is convenient to use Pi + Ci = I to switch 
between Pi and Ci, depending on which operation is easier. To check the equivalence between two visually 
different expressions, in appendix A we introduce a simple trick to solve all constraints and get a set of 
independent kinematic variables. This trick can uniquely fix the form of an expression no matter by which 
means it is obtained (it is better to use this trick with a computer algebra program). 


2.2 Deformation generator algebra 

Now we begin to explore the commutativity of C’s, which can be decomposed into the commutativity at 
integrand level and at integral level. The former is encoded in two successive deformations, and the latter 
is encoded in two successive contour integrals, which will use Laurent expansion m. w = 1/z. Before this, 
we need to first study the BCFW deformation generators and their algebra. 

^The term in the denominator can only originate from spurious pole {i\K\j\, where K contains at least two external 
momenta other than i,j, when it is deformed by {i\j]. All other physical poles can at most contribute terms linear in z under 
one BCFW deformation. 
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Let’s define the BCFW deformation generator with respect to {i\j] as 


D 


(*|j] 




d 


+ Af 


A 


then the familiar BCFW deformation becomes 


( 2 . 12 ) 


exp R{\i, \j) = R{\i - zXj, Xj + zXi). (2-13) 

Although by default the spinorial partial derivatives treat all spinors as independent, we must also impose 
the momentum conservation constraint on real amplitudes. Without doubt, this constraint will affect the 
independence of spinorial partial derivatives, but it will not affect the commutator algebra of . Below 
we will provide a simple argument. 

Note that any automatically annihilates the sum of all external momenta, i.e., 

^ -C'(i|j](AiAj + XjXj) = 0, (2-14) 

so we claim that momentum conservation is a trivial constraint. To get some intuition, one can consider 
a spherical surface, for which any rotation generator, say Lxy, annihilates the constraint 

= r^. (2-15) 


To parameterize one of the spherical symmetries explicitly, we can define an angle 9xy via 

d d d 

^dy ^dx ~ d9xy' 


(2.16) 


while X, y are no longer independent on the sphere, 9xy can be arbitrary, as this degree of freedom moves 
a given point around on a subset of the spherical surface. From this viewpoint, the commutator algebra 
of Lxy, Lyz, Lzx is obviously unaltered. More profoundly, it is these rotation generators that fully generate 
the spherical surface. Given a particular point in M^, rotation generators move it around to sweep over the 
entire surface of a fixed distance from the origin. 

This picture can be exactly generalized to the case of BCFW deformation generators. We can define 
an ‘angle’ in a complex spinorial sense for each deformation, via 




d 

d9^i\j] ’ 


(2.17) 


then parameterizes one of the symmetries that preserve momentum conservation, and hence momen¬ 
tum conservation will not alter the commutator algebra of at all. But instead, this constraint is fully 
generated by 2 = n{n — 1) BCFW deformation generators. Given a particular point in namely the 

complex spinorial space (Aj,Aj), BCFW deformation generators move it around to sweep over the entire 
codimension-4 surface of a hxed sum of external momenta. And physically, this sum is zero. 
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Since the commutator algebra is unaltered, we are free to treat all spinors as independent to derive 
the commutation relations. Imagine the 0th step of deformation is then the 1st step can be one of 
the four types as named below: 


{k\l] = independent, 

(i|/] or {k\j] = straight descendent, 

(2.18) 

(^|i] or {j\k\ = skew descendent, 

{j\i\ = cross descendent. 

The generators of first two types commute with that of (i|j], i-e., 

[-D(j|j], =0, = \DD{k\j^ =0- (2-19) 


For the last two types, 


2hi = —A 




* 9A“ 


dX 


( 2 . 20 ) 


where hi is the helicity operator with respect to the i-th particle, for a function covariant under the little 
group (an amplitude does have this scaling property). 

For the skew descendent case, using the Baker-Campbell-Hausdorff formula 

exp X exp F = exp(X + Y) exp Q [X, F]^ , for [X, [X, F]] = [F, [X, F]] = 0, (2.21) 

and due to the commutativity of a straight descendent pair, we have 

[exp (2o^(j|j]) ,exp (ziF>(j|fc])] = exp (5;o^(i|j] + ^i^{i|fc]) 2sinh (^ZQZiD(^\k^ . (2.22) 


Hence skew descendent deformations {i\j] and {j\k] commute if annihilates the amplitude, however, 
this is a too stringent condition which often trivializes the deformations being used. Therefore in general, 
skew descendent deformations do not commute. 


2.3 Commutativity at integrand and integral levels 

By applying the BCFW deformation generators, we perform one constant extraction on rational function 
R{Xi, Aj) as (different from contours around finite locations, the contour around infinity is clockwise) 

CR= { —R{z) = f— r( (2.23) 


w 


w. 
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where the change of variable is tc = l/z, so the infinity for z is the zero for w. However, the residue at this 
zero is not a naive one. Recall (2.6) in terms of w, it reads 




{hk + bokw)w‘^'^f^ ^ 
{a2k + aikw + aokw'^)‘^’‘ 


+ C0 + Y1 




(2.24) 


^ \^ZK: I I ) I 

a naive substitution of rc = 0 will cause divergence in the third term above. On the other hand, it is clear 
that after the expansion, the third term actually has no simple pole at u) = 0 , since there is already one 
w in the denominator of the integrand. To remove this divergent term, before the contour integration we 
must Laurent expand R{l/w) around w = 0, i.e., we need to first factor out a divergent factor 1/w'^ with 
7 > 1, leaving a finite fraction at tc = 0, then Taylor expand it around w = 0. A simple example is 

-1 

, (2.25) 


bo + biz + 622 ^ 
ao + 012 : 


1 62 + biw + bow^ 


62 


w 


ai + aow 


aiw 


1 I I ^0 2 

1 + —w + —w 

62 02 


1 I “0 
1 H- w 


ai 


then we can Taylor expand the finite fraction around rc = 0, and the contour integral will only pick up the 
constant part in this expression. Similarly, performing two successive constant extractions gives 


with 


CiCoR = 


1 1 


dwi f dwo 


Wi Wo 


0 Jo Wo 


1 


R{—,— 

Wi Wo 


R I —, — = exp —Di exp 


Wl 


—Do 

Wo 


R. 


(2.26) 


(2.27) 


For independent and straight descendent cases, [Di,Do] = 0, so the order of deformations is irrelevant, 
and hence the commutativity of these two types holds at integrand level^. 

However, before performing the integral, to double Laurent expand a fraction is a bit tricky. First, to 
properly factor out the overall factor in terms of wo and wi, we need to ensure that in 

P{wo,wi) 


R( —1 = 

^Wi Wo, 


1 


(2.28) 


wl°wf UiQi(.wo,wi)' 

both P{wo,wi) and Qi{wo,wi) are irreducible polynomials, i.e., there is no common factor in each 

Pi or Qi. Then, we find the expansions of a fraction in opposite orders are different, when any of the Qj’s 
does not contain a constant term. For example, take g{wi,wo) = l/(rci +rco) and expand it around wq = 0 


and for the reverse order, 


i around wi 

= 0 ), 

we 

have 



g{wi,wo) = 

— ( 

/ 

Wo 

< \ 

2 • • • 1 5 

(2.29) 


K 

Wl 

tiif / 

g{wi,wo) = 

— ( 

/ 

Wl 

■ — + 

i 

2 • • • 1 5 

(2.30) 

Wo \ 


Wo 

Wo J 


"‘It is possible that two constant extractions commnte, even if they do not commute at integrand level, but we will not 
consider this trivial case here. 
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hence they are clearly different. In general, if any of the Qj’s happens to satisfy Qi(0,0) = 0, the double 
expansion depends on the order. Conversely, if all Qi’s obey Qi{0, 0 ) / 0 , the order of double expansion is 
irrelevant. Since the contour integral picks up the constant part in the expansion only, the commutativity 
of the denominator expansion is equivalent to the commutativity at integral level. 

In practice, since it is clear that a detectable pole of merely one of the two successive deformations 
always contains a constant term after factoring out a proper factor, we should only focus on the overlap 
of two sets of detectable poles. But it’s impractical to trace each term at each step, for seeing whether a 
constant term exists. 

Combining everything above, one reaches the conclusion; the commutativity of C operators holds for 
independent and straight descendent deformations, provided that the condition Qj(0,0) 7 ^ 0 is satisfied. 
These two types may be considered as ‘good’, as they may enjoy the orthogonal property (2.11). However, 
merely the good types are not sufficient to capture all physical pole terms, as will be explained in the end 
of appendix B. Hence we will not proceed further in this direction, but return to seek for the condition of 
Cn - ■ ■ Cl Co = 0 , such that the (n + 1 ) steps can fully capture the amplitude. 

Nevertheless, this investigation gives a crucial hint for the subsequent analysis: In the expansion, the 
coefficient of after a deformation becomes a new pole of the corresponding boundary term, with order 
one or higher. To be concrete, consider the example below 


11 w 

a a + hz (6 + aw) 


w 


o \ -1 

1 + -W 


w 

~b 


a a 9 
I - -w + - 71 C - 
b 



(2.31) 


where b is the only source of poles afterwards. 


3. Systematic Algorithm of Finite, Definite Steps 

In this section, we propose the systematic process to capture the amplitude after finite, definite steps of 
BCFW constant extractions. The condition of correctly completing the calculation is simply Cn - ■ ■ CiCq = 
0. To achieve this, a form of all poles in the final boundary term is obtained, after a sequence of constant 
extractions which is called the ‘pole concentration’. This sequence is designed for covering all situations so 
there is no extra restriction such as color order, and how to optimize it case by case is set aside temporarily. 
Having the final form of poles, merely using the information of mass dimension and helicities is sufficient 
to judge whether the final boundary term vanishes. 

3.1 Pole concentration 

Now we use pole concentration to capture all poles regardless of whether they are physical or spurious by 
applying the logic of (2.31). Also, each time we perform a BCFW constant extraction on the amplitude, 
at least one of its physical poles will be filtered out, and consequently each corresponding boundary term 
will contain at least one pole mutated from the original physical poles. 




For example, consider denominator (12)(23) (the numerator is neglected for our purpose), under con¬ 
stant extraction (1|3], 

(12)(23) ^ ((12) -z(32))(23) ^ (^’ 

where pole ( 12 ) has been replaced by (32). Crucially, under a next constant extraction, pole (23)^ is either 
unchanged or replaced by another pole of the same order. This means once two poles are stacked, they 
are stacked forever. The same logic also works for anti-holomorphic poles. For a multi-particle pole, we 
first need to turn it into a product of holomorphic and anti-holomorphic poles, with a proper choice of 
deformation. As one example, under constant extraction (1|4], 


1 1 _ 1 

^ Pf23 + ^(4|2 + 3|1] ^ (4|2 + 3|1]’ 

next, under constant extraction (2|5], 

1 1 _ 1 
(4|2 + 3|1] ^ (4|2 + 3|1] -z(45)[21] ^ (45)[21]’ 


(3.2) 


(3.3) 


then we are again left with two-particle poles. 

In general, one can first turn a multi-particle pole into (iijPjji]) where P includes either ii or 
ji, note that pi-^ or pj^ in P is already filtered out by (iij or jjij- Next, one can split {ii\P\ji] by using 
{k\j 2 ] or (^ 2 ^], where P includes j 2 or 12 but not k or 1. This way turns the pole into a product of one 
holomorphic and one anti-holomorphic pole. Then for two-particle poles, once they are stacked, they must 
mutate as a whole afterwards. After finite steps, all poles can be encapsulated in only one holomorphic 
and one anti-holomorphic pole, with orders larger than one in general. 

In appendix B, one sequence of BCFW constant extractions is presented to turn all poles of the final 
boundary term into a common denominator, of the expression given by^ 


(polynomial) 

{ili2)^[hu]^ 


X (remaining factor). 


(3.4) 


where are four different arbitrary particle labels. The remaining factor is a rational function, 

which is dimensionless and helicity-neutral, see (3.8) for example. Note that we have not reduced the 
denominator against the numerator. At the first glance, the reason to get this final denominator is that 
we can use one more deformation, say to get the maximal large 2 suppression, since all poles after 

concentration are vulnerable to it. But in fact, there is a less obvious argument for eliminating the final 
boundary term without introducing one more step, as will be given later. 

Here, m or m gets contribution from physical holomorphic or anti-holomorphic poles, and both of 
them get contributions from physical multi-particle poles. In general, m and m need not be equal, since 
not all possible poles are physical for a particular amplitude. To see the range of m,rn, we will analyze 
all possible physical poles for various n’s. When n = 4, only a half of all two-particle poles can appear 

®The choice of sequence is not unique, and how to optimize it to shorten the steps is a very valuable future problem. 
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in the amplitude, since they are doubly duplicated by momentum conservation. When n = 5, there are 
only two-particle poles, as three-particle poles are equivalent to them by momentum conservation. When 
re > 6, multi-particle poles arise. Their particle numbers range from 3 to (re — 3), to avoid duplications 
of two-particle poles by momentum conservation. To further avoid duplications of themselves, one can 
fix the pole momentum by demanding one pivot particle to be always included, and then the number of 
multi-particle poles is reduced by one half. 

According to the counting above, the maxima of m,m are 

w-max = Tnmax = + • • • + C*” = 2^ ^ — (n + 1), (3-5) 


which nicely covers the special cases of re = 4 and re = 5. 

However, there is a little subtlety in (3.4); For a given amplitude, while m,m can be easily read off 
by analyzing all of its non-vanishing factorization limits, the final boundary term in general contains not 
only poles (^ii 2 )™'[* 3 U]”*) but also the same poles of higher orders from the dimensionless, helicity-neutral 
remaining factor. This also occurs in each intermediate step for each corresponding intermediate boundary 
term. A simple example is the MHV amplitude A(l“, 2“, 3'*', 4^), given by 


,4 - ( 12 )^ 

(23)(34) (41)’ 


(3.6) 


and deformation (1|3] turns it into (recall that z = 1/w) 

((12)-z( 32))3 ^ 1 (-(32) + (12)re;)3 ^ 1 (32)^ / _ ^ ^ 

(23)(34)((41) - z(43)) t(;2(23)(34) -(43)-F (41)'u; t(;2(23)(34) (43) V (32) / V (43) / ’ 

(3.7) 

note that pole (41) is turned into (43), but its order can be larger than one. Explicitly, the corresponding 
boundary term is 


C'(l|3]^ — 


1 

(23) (34) (43) 


7 (41)2(32)3 (41)(12)(32)2 

V (43)2 (43) 


-F 3(12)2(32) 


(12)^(32) 7 (41)2(32)2 _ , (41)(32) \ 

(23)(34)2 1^(43)2(12)2 '"(43)(12) ^ ) ’ 


(3.8) 


where the term in parentheses is the remaining factor of (3.4). The advantage of packing up many pole 
terms into a dimensionless helicity-neutral factor is that, if we can show this representative factor cannot 
exist, the full expression including terms with higher-order poles must also be forbidden. 

One digressive comment is that, so far we have found BCFW deformations to be the only type which 
admits a feasible pole concentration. A counterexample is, there is no straightforward pole concentration 
for Risager deformations [20]. We will not further explain the claim here but it is not hard to confirm it. 
This is another specialty of BCFW deformations, in addition to that BCFW deformations automatically 
preserve (or generate) the momentum conservation constraint. 
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3.2 Kinematic mass dimension 

To prepare for the later analysis, we will study the general information of an amplitude; mass dimension 
and helicities, with which the applicable range of multi-step BCFW recursion relations can be clarified. 

First, for QFTs in 4-dimension, the mass dimension of an n-particle amplitude is (4 — n). We can use 
the LSZ reduction formula to prove this. Schematically, an n-particle amplitude A is defined via 

d^xe^P^eA^ (4>i...<f>„) = A (3.9) 

where (<i>i... is the n-point function, e and A are the wave-function and kinematic operator for each 
field 4>. For a bosonic field, the mass dimensions of e, A and <i> are 0, 2 and 1 respectively, for a fermionic 
field, the mass dimensions of e, A and are 1/2, 1 and 3/2 respectively. Hence the mass dimension of 

J d^xe^P^eA^ (3.10) 

is —1. There are n such pieces, plus the momentum conserving delta function, the mass dimension of A is 
clearly (4 — n). 

One special bosonic field is the graviton. By the perturbative definition g^i, = rj^y -|- it should be 
dimensionless. To treat it as ordinary bosonic fields, we need to redefine it via 


9^iv 


T ^ hfj^y, 


(3.11) 


such that hf^y carries mass dimension 1 and k carries —1. Choosing k to be VSttG, the free field part of 
Einstein-Hilbert Lagrangian takes an analogous form as those of ordinary bosonic fields. Consequently, k 
becomes the coupling constant of gravity. 

One can also rediscover k via the on-shell method. For gravity, three-particle amplitudes including at 
least one graviton are 


^(l_/i,2+h^3-2) ^ ^^^2)-2^23)"^^+^(31)+^'"+^, 

^(l_/i, 2+h^ 3+2^ ^ ^ j^2]-2[23]+2/i+2[3i]-2ft-h2^ 


(3.12) 


where h = 0,1/2,1,2, for all ‘realistic’ theories. No matter which value h takes, k always carries mass 
dimension —1, since the mass dimension of three-particle amplitudes is 1. 

In general, we can reverse the logic above and define the ‘kinematic mass dimension’ of an n-particle 
amplitude as 

D = A - n-'^{Dc)i, (3.13) 

i 

where (Ilc)i is the mass dimension of the coupling constant for each vertex. From now on, we will focus 
on the kinematic part of an amplitude which is a function of spinorial products. 

For the Standard Model®, all coupling constants are dimensionless so H = 4 —n, which is a non-positive 
number for n > 4. For gravitational interactions Dc = — 1 and we will show that D = 2. 


“Strictly speaking, we mean all Standard-Model-type interactions of massless particles, and only dimensionless couplings 


are involved. 
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When D < 0, there is at least one irreducible denominator of the amplitude, which means a pole to be 
detected by BCFW deformations. Conversely, when D > 0, the amplitude may admit some invulnerable 
terms to BCFW deformations which include polynomials and ‘pseudo polynomials’. The classification of 
these objects can be found in appendix C. 


3.3 The master formula 


By using mass dimension and helicities, we now derive the master formula for the subsequent discussions. 
After pole concentration, the final boundary term schematically reads^ 




(3.14) 


(12)m[34]-i^i' 

where we have temporarily taken ii, 2 , 3,4 = 1,2,3,4. The reason to use un-contracted spinors is that, it is 
more compact to capture the helicity information, and it can save the Schouten identity manipulations, as 
one can freely recombine them to get the desired spinorial products. Of course, the cost is that one needs 
to rule out all those illegitimate combinations. This treatment is similar to the methods used in [21, 22]. 
Now the helicity configuration enforces that 


—2hi + m = ai — /3i, 


- 2 /i 2 + m = 02 - /32, 

2/i3 + m =/33 - as, (3.15) 


2/i4 + m = (34 — 04, 

2hi = I3i- Oi, (i = 5,... ,n) 

where m,m are known for a particular amplitude. Note that there are 2n variables, with only n helicity 
constraints. We will fully exploit the n remaining degrees of freedom to derive the master formula. 

The kinematic mass dimension of (3.14) is 

1 / n ” \ 

=-(m + m) + - I + I , (3.16) 

^ Vi=l i=l J 

where obviously, ^ o and ^ (3 must be both even to form spinorial products. Also, we have m + m > 1 
with 7T7-,m > 0, and a,/3 > 0. For a legitimate final boundary term, D' equals to D defined in (3.13). 

When D ^ D' under all circumstances, the correct dimension and helicities cannot be satisfied simul¬ 
taneously, and then the final boundary term is eliminated. One direct way to achieve this inconsistency is 
to show is larger than D. First we need to figure out this minimum by eliminating one variable for 
each particle, as there are two variables Oi and (3i to be chosen. For i = 5,... ,n, 

— (aj -|- (3i) = —hi + (3i = hi + Oi, (3-17) 

^In general, the final bonndary term’s numerator is a polynomial but we only focus on one term, as the identical analysis 
applies to all terms. The remaining factor of (3.4) is dropped, since it does not contribute to the master formula. 
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when hi is negative, {—hi + /3i) is guaranteed to be positive, similarly when hi is non-negative, {hi + ai) is 
guaranteed to be non-negative. To manifest the non-negativity of D', our choice is 


^{ai -h Pi) = \hi\ -h mm{ai,Pi). 


(3.18) 


Extending this logic for all particles, yields 

-hmin(aj,/3j)j -h ^ 


D' = — (m-l-m)-l- ^ y^hi 
1 = 1,2 


m 


i=3A 


, m 
hi + — 


-F mm{ai,Pi) j -F Ed'" j| -hmin(aj,/3j)), 

(3.19) 


i=5 


which is the master formula, and explicitly, 


'^{\hi\ + min(ai,/3i)) = '^{-hi + Pi) + '^{hi + a*), (3.20) 

i=5 h<Q h>0 


which separates the sum into two parts according to the helicities. The final boundary term (3.14) now 
reads {pi = |z)[i| is a helicity-neutral momentum with additional mass dimension 1 ) 


_rr 

(12)™[34]™ 




{i 


—2/ij+m 


Pi 


n 

i=3A 


i\2hi+m a, 

0 Pi 


{i 


-2hi-m 0'. 


n ('I ■“*»'?* ni'i 


2hi 


^Pi 


h<0 


h>0 


(3.21) 

where / means one of two candidate expressions is chosen to manifest the non-negativity of D', as this 
choice also manifests the ‘extra neutral momenta’, in addition to the ‘net spinors’ that carry the helicity 
information. While the latter content is mandatory, the former is optional since it is brought in to fill the 
extra capacity of mass dimension. There is no unique choice of picking these extra a’s and P’s as long as 
the total dimension is correct. 

For i = 5,..., n, \hi\ is trivially non-negative. For hi, /i 2 , /13, / 14 , careful analysis is needed as m, m are 
involved. Rewrite (3.19) as 


n n 

D' = -{m + m) + ri234 + ^ \hi\ + '^Tain{ai, Pi), 

i=l i=\ 


where 


^1234 = ( 

, m 

hi - 

2 

- I^il) + X] ( 

, m 

hi + - 

i=l,2 


2=3,4 ^ 




since m,m and ^ \h\ are fixed, we only need to manipulate ri 234 - It’s easy to check that 


(3.22) 


(3.23) 


m 


hi - 

2 


-\hi\ 


m 


hi + — 


-\hi 


m/2, /ij < 0 

< m /2 — 2hi, 0 < /ij < m /2 

—m/2, m/2 < hi 

—m/2, hi < —m/2 

< m /2 -|- 2hi, —m/2 < hi < 0 

^ m/2, 0 < /ij 


(3.24) 
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i=l,2,3,4 


Figure 1: Sum of the distances of four line segments. It is obviously no less than (m + m). 


to maximize ri 234 , one must take to be two of the smallest helicities, and to be two of the 

largest helicities in the process®. 

On the other hand, even if one chooses /ii,/i 2 ,arbitrarily among all hj’s, is no less than 

zero after a similar relabeling. To see this, let’s rewrite (3.19) as 


n n 

D' = -{m + m) + r (234 + ^ \hi\ + ^ min(ai,/3i), 

2 = 5 2=1 


and focus on the quantity 


T' = 
^1234 — 


\hi 



i=3,4 



(3.25) 


(3.26) 


by this definition T (234 has a simple geometric meaning: It is the sum of the distances of four line segments 
stretching from /ii,/i 2 to line h = m/2, and from h^,hi to h = —m/2, as shown in Figure 1. It’s easy 
to find that its minimum is (m + m), when four points are on one horizontal line and the line is within 
the region between h = m/2 and h = —m/2. When this horizontal line moves outside the region, ^(234 
increases by 4 x (distance above or below the region). When this line is not horizontal, one can see that 
^1234 i® always larger than (m + m) after a relabeling such that hi, /12 < h^, h^. Since T [234 is no less than 
{m + fn), T*(njn is always non-negative. 


®Now we should strictly use 0 , 12 , * 3 , *4 instead of 1,2, 3, 4 to admit a possible relabeling, when we implement the desired 
arrangement of pole concentration. 
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Also note that 






\hi\ = Oy/ {-2hi) j (-m), 
\hi\ = {-m)j {2hi) j 0 , 


(3.27) 


which always take integer values. According to (3.22), if |/i| is fractional, D' must be fractional. 


3.4 Inconsistency elimination 

We are prepared to show that all massless tree amplitudes except those admitting (pseudo) polynomials of 
given mass dimension and helicities, shall be fully determined by multi-step BCFW recursion relations. In 
the following analysis, we pretend not to know any knowledge of QFT in the Lagrangian paradigm besides 
mass dimension and helicities. These are the only data needed to construct three-particle amplitudes, and 
recursions extend them to all higher-point ones. 

First note that, after pole concentration, there is no need to find any further deformation to kill the 
final boundary term, because another choice will not change its pole form, so one can always rearrange 
the entire sequence of pole concentration to reach the desired relabeling. Hence it must be the last step if 
chosen properly, and the direct way to use is the ‘inconsistency elimination’. 

Within its own framework, if inconsistency elimination can exclude the final boundary term, the new 
algorithm should be a completely independent approach for calculating a particular amplitude. Otherwise, 
terms that survive it need to be identified, similar to (pseudo) polynomials. In fact, we do discover a new 
type of object called the ‘saturated fraction’ in this way. 

From the master formula (3.19), it is already known that > 0. Therefore if H < 0, H 

always holds. This tells us the nontrivial cases are of H > 0. When D > 710^^, we have the inconsistency 
criteria below to further eliminate the final boundary term: 


(1) Fractional Dimension (FD): If D' = fractional. This arises when ^ \h\ = fractional, but Lorentz 
invariance demands the dimension of (3.14) to be an integer. This also implies that fermions must appear 
in pairs to be consistent. 


(2) Pair Mismatch (PM): If ^ a = odd and X]/3 = odd, when FD is excluded already. In this case, 
(3.14) cannot be written as a fraction in terms of Lorentz invariant spinorial products, even though the 
dimension of (3.14) is an integer. 

(3) Spinor Excess (SE): If there exists an f, such that ctj > or Pi > case, 

spinorial contraction will force (3.14) to vanish, even though a = even and /? = even. 

Altogether, there are four layers of inconsistency criteria: 

(0) D < 
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(1) If Z) > consider FD. 

(2) If and FD is excluded, consider PM. 

(3) If D > and both FD and PM are excluded, consider SE. 

It’s obvious that these inconsistency criteria only mention general properties of field theories. Hence 
inconsistency elimination is theory independent, while in practice knowing some theory dependent proper¬ 
ties would help simplify discussions case by case. If the final boundary term can survive all four criteria, it 
must admit ‘saturated fractions’ (SF). Note that we have already set aside (pseudo) polynomials, because 
they can be identified without using the master formula (3.19). Altogether, there are three types of objects 
invulnerable to BCFW deformations: 

(1) Polynomials, such as (12). 

(2) Pseudo polynomials of n = 4, such as [34]/(12). When (12) —>■ 0, we must also have [34] —>■ 0 since 
we are using a BCFW deformation. Then the ratio [34]/(12) is finite like a polynomial. 

(3) Saturated fractions of n > 5, such as [34] [56]/(12) of n = 6. When (12) —>• 0, the fraction becomes 
divergent, so it is different from pseudo polynomials. But somehow similarly, any BCFW deformation of 
n = 6 fails to render such a fraction vanish under the large z limit. 

Among these three objects, a polynomial is completely inert to BCFW constant extractions, in fact 
it is invulnerable to any type of deformation in on-shell methods. A pseudo polynomial is also completely 
inert to BCFW constant extractions, while this requires momentum conservation. A saturated fraction is 
form-inert to BCFW constant extractions®, but with particle labels rearranged. The latter two objects are 
vulnerable to other types of deformations, such as Risager deformations [20]. The detailed exploration of 
all these three types is presented in appendix C. 

So far, we have witnessed how the systematic process of multi-step BCFW recursion relations can be 
arranged to determine a particular amplitude. In summary, there are four steps: 

(1) Analyze all non-vanishing factorization limits to determine the amplitude’s common denominator, 
which is a product of all physical poles. This stage can be done almost purely diagrammatically. 

(2) Figure out the amplitude’s kinematic mass dimension, then combine this with its helicity configu¬ 
ration to identify all possible (pseudo) polynomials. If none of them arises, we assume the amplitude can 
be fully determined and proceed to the next step. 

(3) Choose four particle labels (two of the smallest helicities and two of the largest ones) to determine 
the denominator’s form of the final boundary term, and arrange a sequence of BCFW constant extractions 
to proceed pole concentration. This sequence must be able to capture all physical poles, such that each 
of them contributes to the final denominator via powers m,m. Ensuring this, the sequence should be as 
concise as possible. Such an optimization is a very valuable future problem. 

(4) Use all four inconsistency criteria layer by layer to eliminate the final boundary term. If it fails, 
identify all possible saturated fractions. Then discuss whether these saturated fractions are legitimate, if 

®Here we mean a pure saturated fraction. A mixed saturated fraction is a pure saturated fraction times a polynomial. Its 
transform under a constant extraction will be demonstrated in appendix C. 
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not, clarify the argument to rule them out as described in the end of appendix C. This delicate treatment 
to remove all dependence on spurious poles is another valuable future problem. 


4. Applications in Standard Model plus Gravity 


Knowing the general guide of multi-step BCFW recursion relations, naturally we would like to see how it 
applies to specific theories. As familiar examples, realistic theories, i.e., the Standard Model plus gravity^® 
are being considered. For simplicity, we will assume that all particles are massless. 

For reader’s convenience, we rewrite the master formula below 


D' = —{m + m) + ^ \hi 
i=l,2 


m 

y 


+ E 

2=3,4 



n n 

+ X] 1^*1 + X] A ); 

2=5 2=1 


(4.1) 


recall that one should take to be two of the smallest helicities, and to be two of the largest 

helicities in the process. 


4.1 Tw^o separated sectors 

Let’s first consider Standard Model and gravity separately. For the Standard Model, the previous section 
gives D = A — n. On the other hand, since is no less than zero, when H < — 1, the final boundary term 
must be eliminated. This directly tells that all Standard Model amplitudes of n > 5 are solvable, leaving 
amplitudes of n = 4. The n = 4 case admits (pseudo) polynomials 1, ((34)/[12])^^ and ((34)/[12])^^, but 
the last one will be excluded. As it is proved in the next subsection, any amplitude’s helicity configuration 
in gauge theory must be between MHV and anti-MHV. 

For pure gravity, assuming that one of the Feynman diagrams of an n-particle amplitude contains Vm 
m-point vertices and p internal propagators, it is clear that 


'^mvm-2p = n, '^Vm=P + l, '^{m - 2)vm = n - 2, (4.2) 

mm m 

and each m-point vertex brings in (m — 2) k’s, hence from (3.13) we have 


D = A — n — (—1) '^{m — 2)vm = 4 — n -|- (re — 2) = 2, 

m 


(4.3) 


now we compare D with D'. When re > 6, = 4 > 2 so this is completely solvable. When re = 5, from 

(4.1), and in the most conservative situation —rre/2 < —2 < 2 < m/2, the all-plus helicity configuration 
(similar for all-minus) admits the saturated fraction 


[34] 2 [35] 2 [45] 2 


(4.4) 


^°Any massless theory can be analyzed analogonsly. However, as the amplitude’s kinematic mass dimension goes up, more 
types of (pseudo) polynomials and saturated fractions may arise, and one needs to identify them carefully. 
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When n = 4, the all-plus helicity configuration admits pseudo polynomials such as 

If (13)<24). (4.5) 

where x, y are unspecified, while the all-but-one-minus case gives = 4 > 2 already. 

However, similar to gauge theory, any amplitude’s helicity configuration in pure gravity must also be 
between MHV and anti-MHV. And for the MHV configuration, (4.1) gives = 8 > 2. Therefore pure 
gravity is completely solvable. 

In general, for Standard Model plus gravity we have D <2. Note that an amplitude which contains 
gravitational vertices only always obeys D = 2, regardless of how many or what kinds of external legs it 
owns, as later shown by (4.10). This specialty implies that one can arbitrarily attach more particles to a 
known amplitude without changing its mass dimension, via gravitational interactions. 

4.2 MHV configuration of gauge theory and gravity 

For the usual gauge (Yang-Mills) theory and gravity, based on the knowledge of multi-step BCFW recursion 
relations, one can prove that the helicity configuration of any non-vanishing amplitude must be between 
MHV and anti-MHV, i.e., there is no all-plus or all-but-one-minus configuration. 

From Lorentz invariance and little group scaling, the three-particle amplitudes for gauge theory and 
gravity are known to be (s = 1 , 2 ) 

A(l-2-3«) ^(1-2-3-) (4.6) 

while the T^-type (s = 1 ) or the i?^-type (s = 2 ) three-particle amplitudes are 

A(l-^ 2-* 3-*) = gs, _((12)(23)(31))* , A(l+^ 2+^ 3+^) = gs,+++ ([12] [23][31])^ . (4.7) 

For gauge theory the coupling constant is dimensionless and for gravity the coupling constant carries mass 

dimension —1, while mass dimensions of gs, _and 55 ^+++ are —2 for s = 1 and —5 for s = 2, so the 

all-plus and all-minus three-particle amplitudes are excluded. 

To exclude the all-plus and all-but-one-minus amplitudes, we need to show that these configurations 
have vanishing contributions from either (BCFW detectable) factorization limits, or (BCFW undetectable) 
invulnerable objects which include (pseudo) polynomials and saturated fractions. Recall the key identity 
of multi-step BCFW recursion relations 

I = Pn CnPn-l + . . . + CnCn-1 ' ' ' <7*2^1 -|- CnCn-1 ' ' ' C 2 C 1 PQ + CnCn-1 ' ' ' C'2C'iC'o, (4.8) 

the terms with Pi's and the term Cn ■ ■ ■ CiCq represent these two parts of contributions respectively, so 
when they both vanish, the corresponding amplitude must vanish. 

For factorization limits, an inductive observation is: An all-plus amplitude can only factorize into one 
lower-point all-plus amplitude and one all-but-one-minus amplitude. An all-but-one-minus amplitude can 
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Figure 2: 4-point gauge vertex is ‘equivalent’ to two connected 3-point vertices. Wavy lines represent gauge bosons. 


only factorize into two lower-point all-but-one-minus amplitndes, or one lower-point all-plus amplitude and 
one MHV amplitude. 

For the all-plus case, we will finally recnrse down to ^(-|- + +), so its factorization limit is zero. This 
also excludes the second possibility of factorization for the all-but-one-minus amplitude, then for the latter 
we will finally recurse down to 74(-|- -|- -|—), which can further factorize into two ^(-|- H—)’s superficially. 
However, when we use a BCFW deformation to calcnlate this part, its contribution is zero. This is due to 
the fact that both of its sub-amplitudes are anti-holomorphic functions, while the non-vanishing BCFW 
constrnction requires one to be holomorphic and the other anti-holomorphic. 

Therefore (BCFW detectable) factorization limits of the all-plus and all-but-one-minus configurations 
are zero, then we consider (pseudo) polynomials and saturated fractions. Note the previous discussion has 
not separated gauge theory and gravity yet, hence it holds for both. 

For gauge theory, the only invulnerable object is the pseudo polynomial of n = 4, as the reader may 
check this in appendix C, given by 

(4.9) 

which is related to A(-|- -|- -|--|-). However, when one nses the holomorphic factorization limit (which is the 
only type of effective deformation to detect pole (12)^, snch as Risager deformation as a familiar example) 
to send Ai —A 2 , this pole is spurious since its order is two, which excludes the pseudo polynomial above. 

For gravity, the only invulnerable objects are saturated fraction (4.4) related to A{+ and 

pseudo polynomial (4.5) related to H(-|- -|- -|-+). Since both of them contain spurious poles, the argument 
above can also excluded them. Now the proof is done. 

We would like to emphasize that here we have used general information only, namely mass dimension, 
helicities and factorization limits, without any further aid such as supersymmetry. 



4.3 Simplified diagrammatic rules 

To simplify the general discussion of Standard Model plus gravity, we introduce the diagrammatic rules 
called ‘stretch and shrink’. The first example is the gauge interaction, as shown in Figure 2. 

Fixing four external gauge bosons, this 4-point vertex can be ‘stretched’ into two connected 3-point 
vertices, without changing the vertex’s mass dimension. This tricky equivalence holds at the level of mass 
dimension and helicities, which are the only information required for inconsistency elimination. 
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Figure 3: Stretch or shift rule of a gravitational vertex. Bold lines represent Standard Model particles, and zigzag 
lines represent gravitons. 



Figure 4: Gravitational component and its trivial case. The unspecified lines between two external lines can be of 
either type of these two specified lines ‘on the boundary’. In this case, they can be either Standard Model particles 
or gravitons, or a mixture of both. 


In other words, we have chosen a representative sub-diagram to encode the same information of mass 
dimension and helicities, and reduce the types of equivalent sub-diagrams in the analysis. Following this 
logic, all higher-point vertices in Standard Model plus gravity can be stretched into a number of connected 
3-point vertices, except the special (j)^ vertex. This simplified rule is notably advantageous in gravitational 
interactions. 

As shown in Figure 3, gravitons can be shifted from any place to any place in a sub-diagram, without 
changing its mass dimension. Physically, this is because gravity is universal, gravitons can emit from any 
part of a system. Mathematically, this is because an m-point gravitational vertex carries coupling constant 
k carries mass dimension —1. This vertex can contain gravitons only or have Standard Model 
lines attached to it. Therefore, the m-point vertex can be stretched into (m —2) connected 3-point vertices, 
with the exception of vertex. 

For convenience we define the ‘gravitational component’, as shown in Figure 4. All vertices within this 
component are gravitational, while its external legs can be either Standard Model particles or gravitons, 
or both. There is one special graviton which will attach to another component. A trivial case is that there 
is no vertex at all, so this special graviton becomes the only component. 

Gravitational components also obey the simplified rules, and for convenience they are usually shrunk 
into one component, as shown in Figure 5. This pack-up can reduce many sub-diagrams of gravitational 
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Figure 5: Gravitational components can be packed into one component. 




Figure 6: All representative vertices in Standard Model plus gravity. Solid lines represent fermions, and dashed 
lines represent scalars. 


components to one sub-diagram. It’s free to attach (or detach) a gravitational component to (or from) 
another component, since this operation will not change its mass dimension. 

Summarizing the simplified diagrammatic rules, we are now left with the representative vertices only, 
as shown in Figure 6. 

4.4 Amplitudes and (pseudo) polynomials of ZD = 0, 1, 2 

In Standard Model plus gravity, the nontrivial cases are ZD = 0,1, 2. By using the representative vertices 
in Figure 6, the following discussion is considerably shorten. 

D = 0 case: First we consider ZD = 0, all possible amplitudes are listed diagrammatically in Figure 7. 
Similar to gravitational components. Standard Model components presented here only contain Standard 
Model vertices. Also note the Standard Model components attached by a single graviton and a nontrivial 
gravitational component, are listed separately for clarity. 

For ZD(jjjj^ = 0, (pseudo) polynomials arise when all helicities are the same. Now the last three diagrams 
in Figure 7 are excluded, since a three-point Standard Model vertex can never have three same helicities. 
The second diagram is also excluded since a graviton has helicity ±2. 
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Figure 7: Amplitudes oi D = 0. 



Figure 8: (Pseudo) polynomials oi D = 0. 


Therefore, (pseudo) polynomials come from the hrst and third diagrams, as listed in Figure 8. The 
first three diagrams in Figure 8 admit polynomial 1, while the fourth diagram admits pseudo polynomial 
([34]/(12))^^, as the Yukawa interaction permits the two fermions of its vertex to have the same helicities. 
The fourth diagram cannot have a gravitational component attached to it while maintaining the four same 
helicities, since fermions and gauge bosons must appear in pairs of opposite helicities when coupling with 
gravitons. Finally, since four gauge bosons must take the MHV conhguration, ([34]/(12))^^ is excluded. 

All amplitudes oi D = Q diagrams other than those in Figure 8 are solvable. We find it convenient to 
attach a gravitational component to a known diagram to simplify the discussion. This one-line attachment 
is equivalent to maximally separating gravitational and Standard Model vertices into two sub-diagrams, 
when building the simplest representative diagram. 

D = 1 case: Continuing in the same fashion for D = 1, all possible amplitudes are listed in Figure 9. 
Here, the 3-point Standard Model vertex can be one of the following four types: (a) 3-gauge interaction 
(±1,-)-l, —1); (b) gauge-fermion-fermion interaction (±1,-|-l/2, —1/2); (c) gauge-scalar-scalar interaction 
(±1,0,0); (d) scalar-fermion-fermion (Yukawa) interaction (0, ±1/2,—1/2). 

For the left diagram in Figure 9, when vertices of type (a), (b), (c) and (d) are attached by a gravi- 
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Figure 9: Amplitudes of D = 1. 




Figure 10: Polynomials of D = 1. 


ton, its four helicities are (±1,+1, —1, ±2), (±1,+1/2,—1/2, ±2), (±1,0,0, ±2) and (0,+1/2,—1/2, ±2) 
respectively^^. Plugging the data in (4.1), corresponding H+^^’s are 3, 3, 3 and 2, which exclude all four 
cases. In one words, the left diagram is excluded simply due to the single graviton. 

For the right diagram in Figure 9, when vertices of type (a), (b) and (c) are attached by a gravitational 
component, in the most conservative situation this component only contains external scalars, since higher- 
spin particles must appear in pairs of opposite helicities, which will not decrease Its helicities are 

(±1,+1, —1,0,0,...), (±1,+1/2, —1/2,0,0,...) and (±1,0,0,0,0,...) respectively, and ... denotes more 
scalars besides the minimal five. Applying (4.1), corresponding H+^^’s are 3, 2 and 1, which exclude first 
two cases. However, the third case is also excluded even if its is allowed. The argument is that no 
spinorial product can be formed by only |1)^ or |1]^, which is known as the Spinor Excess of inconsistency 
elimination. The only polynomial comes from the vertex of type (d), as given in Figure 10. This polynomial 
is (12) or [12]. 

D = 2 case: The last case is H = 2. The possible amplitude is given in Figure 11, and corresponding 
polynomials are listed in Figure 12. The first one is P^y, where x, y are two unspecified scalars. The second 

^^Here ±1 and ±2 are independent, they do not necessarily have the same sign. 
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Figure 11: Amplitudes of D = 2. 




one is (la:)[x2], with one pair of fermions of opposite helicities. The third one is (12) [34], with two pairs of 
fermions of opposite helicities. One may consider a fourth one, with one pair of gauge bosons of opposite 
helicities, which is allowed since its is 2. But this case is also excluded, as no spinorial product can 
be formed by only |1)^|2]^ or |1]^|2)^. 

Note that these D = 2 polynomials can be of either n = 4, or n > 5 for which all unspecified particles 
are scalars. For n = 4, there are dimensionless pseudo polynomials of the form ([34]/(12))*, which can be 
an additional factor of the polynomials above. This factor will lead to a global shift of all four helicities. 
But incidentally, there is no extra legitimate pseudo polynomial. 

Last but not the least, we need to check the analysis above has covered all possible diagrams, by using 
a compact formula 

5 

D = 2-Y,{si-2), (4.10) 

i=l 

where D is the kinematic mass dimension of a Standard Model ‘skeleton’, as shown in Figure 13. This 
skeleton amplitude is made of S Standard Model components connected by internal gravitons and s* > 3 
is the number of external legs for each component. When Sj = 2, there is no Standard Model vertex, and 
each component reduces to a Standard Model line, which should be excluded by the skeleton’s definition. 
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Figure 13: Standard Model skeleton. In this case, the unspecified lines can be Standard Model particles only. 


The proof of this formula is simple. For S connected components, there are {S — 1) internal gravitons. 
Each internal graviton has two attached points, which bring in two k's. Hence from (3.13), we have 

s 

D = A-Y,Si-{-2){S-l), (4.11) 

i=l 

which is identical to (4.10) after a simple rearrangement. 

Having this skeleton, to build a general amplitude, more gravitational components (either nontrivial 
components or single gravitons) can be attached to it. By applying the simplified diagrammatic rules, they 
can be packed into one single gravitational component. 

Now let’s consider D = 0, we have S' = 1, si = 4 and S = 2, si = S 2 = 3. These two cases correspond 
to the hrst and fourth diagrams in Figure 7. When attached by gravitational components, they become the 
rest four diagrams. For D = 1, we only have S = 1, si = 3. This diagram cannot stand alone, since there 
is no physical massless 3-particle amplitude. Hence it must be accompanied by gravitational components, 
which gives the two diagrams in Figure 9. For D = 2, the amplitude contains only gravitational vertices, 
which corresponds to the diagram in Figure 11. Therefore all possible diagrams have been covered. 

5. Discussions 

In this work, we are mainly concerned with the workability of multi-step BCFW recursion relations. The 
key techniques of this approach are pole concentration and inconsistency elimination. Its applicable range 


- 25 - 






has been clarified and we find three types of objects invulnerable to BCFW deformations: polynomials, 
pseudo polynomials and saturated fractions. While the last two objects can be determined by other types 
of deformations, how to deal with the first one is probably beyond usual on-shell methods and it may lead 
to important generalizations of the present approach. Moreover, when saturated fractions arise, we need 
to discuss whether they are legitimate, if not, finding an argument to rule them out is a valuable topic. 
Again, we would like to emphasize that this systematic algorithm uses general properties of field theories 
only, such as Lorentz invariance, locality and unitarity. The major information we use are mass dimension 
and helicities. 

Ensuring its workability, we try to further improve the efficiency of multi-step BCFW recursion rela¬ 
tions by taking two sophisticated aspects into account, as listed below: 

(a) Knowing the (final or intermediate) boundary term’s schematic form (in terms of un-contracted 
spinors), it is also very natural to seek for a deformation which renders the boundary term vanish under 
the large z limit, other than employing inconsistency elimination. In practice, one can consider both ways 
at each step to shorten the sequence of pole concentration. More profoundly, inconsistency elimination is 
only an argument afterwards, as any boundary term that vanishes must be killed by a good deformation 
with respect to that step. 

(b) In practice, it is often evident that there is no need to reach the final denominator (3.4). Merely a 
particular intermediate form is sufficient to complete the calculation correctly. This is due to the fact that 
pole concentration is for eliminating as many neutral momenta as possible in the denominator, and hence 
increasing the ‘net spinors’ in the numerator while fixing the helicities. 

Let’s illustrate these two aspects through three simple examples. The first one is the MHV amplitude 
A(l-,2+,3-,4+,... ,??-■*■), where 1 and 3 are non-adjacent. Assuming all lower-point MHV amplitudes 
are known, non-vanishing factorization limits give the common denominator 


(12)(23)...(n-l,n)(nl)’ 

on the other hand, the kinematic dimension of this amplitude is (4 — n), as the gauge coupling constant 
is dimensionless. The correct helicities require its numerator to be |1)^|3)^ schematically, which uniquely 
fixes the amplitude as 

( 13 )^ ,5 2 , 

(12){23) ...(n-l,n){nl)’ ' ' 

from this example, we see the schematic form is a simple but powerful tool. 

The second example is amplitude 2"’'^, 3“^, 4'’'^, 5“^) in Einstein-Maxwell theory. Non-vanishing 

factorization limits give all physical poles as [12] [32] [14] [34] [15] [25] [35] [45]. Since its kinematic dimension 
is 2, the correct helicities require its schematic form to be 


[1| |2|^|3| [4|^xn*h>[»l 
|12|[32][14]|34][15][25][35][45]’ 

where ]•)[•]’s are unspecified neutral momenta. At this stage, one can already check that 


(115] (215] (315] (415] 


(5.4) 
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are good deformations since each of them induces in the denominator, while the numerator can at most 
contain z'^ to avoid Spinor Excess, i.e., j*) in the numerator can at most contain two identical spinors 

to form non-vanishing spinorial products, which restricts n^[*l ™ the same way. For this case, to find 
a deformation with maximal large z suppression is clearly more straightforward than using inconsistency 
elimination after a number of constant extractions. 

The third example is even more interesting. Consider amplitude 3“^, 4+^, 5“^, 6"*“^) again 

in Einstein-Maxwell theory. Factorization limits, mass dimension and helicities together fix its schematic 
form to be 


(lp[2P(3|'"[4p(5r-!|6|'"xn“K>l' 


p2 p2 p2 p2 p2 p2 p2 p2 p2 ^ 

^ 12 ^ 14 ^ 16 ^ 32 -^ 34 ^ 36 ^ 52 ^ 54^56 ^ 


p2 p2 p2 p2 p2 p2 p2 p2 p2 ’ 

^ 132 -^ 134 -‘^ 136 -^ 152 -‘^ 154 -‘^ 156 -‘^^ 352 -^ 354^356 


and we will show that two successive deformations, namely 


(5.5) 


(3|1]^(5|1] (5.6) 

can already capture the full amplitude. First, after constant extraction (3|1], similar trick of pole concen¬ 
tration gives its schematic form 

_ (l|^l2p[3|^0[4p(5p[6pxn^^|.)[.| _ 

^52^l4n'6^i32^f3 A(12)2(14)2(16)2[32] + 2|3]2(1|5 + 4|3]2(1|5 + 6|3]2 ’ ^ ^ 

and constant extraction (5|1] turns it into 

_ (lp°[2p[3|^°[4p[5|^[6pxn^^|.)[.| _ 

(12)3(14)3(16)3[32]2[34]2[36]2[52][54][56](1|3 2|5](1|3 4|5](1|3 6|5](1|5 2|3]2(1|5 -h 4|3]2(1|5 -h 6|3]2 ’ 

(5.8) 

then there is no need to proceed further, because in the numerator Spinor Excess already arises, as !•) 
can never saturate (Ip® to form non-vanishing spinorial products. By this way, two steps can already get 
the correct answer, while a blind pole concentration in general requires 4(6 — 3) = 12 steps. Therefore, it 
is not always necessary to reach the final denominator, when eliminating part of neutral momenta in the 
denominator enforces the numerator to contain sufficient net spinors for triggering Spinor Excess. 

But this is not the end of the story. When proceeding with the calculation 


I — -P(5|l] + C'{5|l]-P(3|l] + C'(5|1] (7(311], (5.9) 

while it is just shown that (7(5|i](7(3|i] = 0, incidentally we also find (7(5|i]P(3|i] = 0. This means (5|1] is 
a good deformation and hence one step is enough. Since particles 3“i and 5“^ are symmetric in the 
helicity configuration (there is no color order), (3|1] is also a good deformation. 

In general, there is a last good deformation corollary: After the n-th step, when I = (known terms)„ -|- 
Cn Co is reached, we can further expand it by one more step as 


I = Pn+i + (7n+i (known terms)^ Cn+iCn ■ ■ ■ Cq 


(5.10) 



assume the (n + l)-th step is the last step, for which Cn+iCn • • • Co = 0, and if 

Cn+i (known terms)„ = 0, (5-11) 

then the (n + l)-th step is a good deformation. This corollary is powerful in practical calculations, since 
unnecessary steps can be saved if we incidentally encounter the condition above. 

Back to the mainline, these two aspects (a) and (b) will be demonstrated more systematically, with 
more examples in our future work, with a possible joint use of the last good deformation corollary. The 
major goal is to improve the efficiency provided the workability is ensured. The exit of this maze is now 
found, and how to shorten the correct route is a complicated yet fascinating problem. Finally, we would 
like to highlight the power of simple analysis by mass dimension and helicities, as these cheap information 
possibly lie in the core of the future study of efficiency. 
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A. Independent Kinematic Variables 


In the calculation of amplitudes, it is common that two visually different expressions in terms of spinorial 
products are in fact equivalent. Although one can use a numerical method to check this equivalence, it is 
still favorable to find an analytic way. 

For an n-particle amplitude, let’s start with all holomorphic spinorial products as listed below; 


( 12 ) 

(13) 

(14) 

(15) 


(23) 

(24) 

(25) 


(34) 

(35) 


(45) 


(A.l) 


(1, n — 1) (2, n — 1) (3, n — 1) (4, n — 1) ... (n — 2, n — 1) 

(In) (2n) (3n) (4n) ... (n —2,n) (n —l,n) 

We can pick |1), |2) as two reference spinors, and for i,j ^ 1,2 all independent Schouten identities can be 
solved via 


(u) = 


1 


( 12 ) 


(li) (2i) 

(l j) (2j) 


(A.2) 


There are (u)’s and C '^_2 Schouten identities, so there are — C '^_2 = 2n — 3 independent {ij)’s. 
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We can repeat this for anti-holomorphic spinorial products, and get (4n — 6) independent kinematic 
variables as listed below: 


(12) 


[12] 


(13) 

(23) 

[13] 

[23] 

(14) 

(24) 

[14] 

[24] 

, n — 1) 

(2,n-l) 

[l,n- 

1] [2,n- 

(In) 

(2n) 

[In] 

[2n] 


But the momentum conservation has not been imposed yet. Adding this constraint, we can solve for [13], 
[14], [23] and [24] for example, via 


(2] j;Pll]=0: 
(1] J;P12]=0: 
(21^P12]=0: 


(13)[13] + (14)[14] = -Sn-(12)[12], 
(23)[13] + (24)[14] = -S2i, 

(13)[23] + (14)[24] = -Si2, 

(23)[23] + (24)[24] = -S 22 - (12) [12], 


(A.4) 


where [i^]- Now there are (4n — 10) independent kinematic variables left. Any expression 

in terms of the remaining variables is uniquely fixed. 


B. Example of Pole Concentration 

Here, we introduce one systematic sequence to turn all poles of the final boundary term into the common 
denominator (ii^ 2 )™[^ 3 * 4 ]™'- la this example, there are four series of BCFW constant extractions, namely 
I, II, III and IV: 


I : 

(2|3] 

(3)4] 

(n — 2 n 

-1] 

(n — In] 

II : 

(l|n- 1] 

(2ln-l] 

(n — 3 n 

-1] 

(n — 2 n — 1] 

III : 

(3|2] 

(413] 

(n — 3 n 

-4] 

(n — 2 n — 3] 

IV : 

(n-311] 

(n-3l2] 

(n — 3 n 

-5] 

1 

e 

CO 

1 


where III and IV and can be copied from I and II, by reducing n to (n — 2) and swapping the holomorphic 
and anti-holomorphic deformed spinors. They have 4(n — 3) steps in total. 

To see how this fully works, we assume that all physical poles occur (restrictions such as color order, 
must be disregarded). For convenience, let’s define sets H = (1,... ,n) and A = [1,... ,n] to denote all 
holomorphic and anti-holomorphic poles of n particles respectively. Note that H and A are not fixed, but 
change at each step of deformation. 
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Also, the set of all multi-particle poles is denoted by M. To fit the analysis of pole concentration, we 
will classify all multi-particle poles according to which particle is absent, in a default order. To be concrete, 
they are listed as 


(^ 2 +^ 23 )^ {P23 + X24)‘^ {P2A + X25)‘^ ... (P2,n-4 +-A2,n-3)^ (-P2,n-3 +-A 2 ,n- 2 )^ (.P 2 ,n -2 +-A 2 ,n-l)^ 

(B.2) 

where Pij = pi + ... + pj, Xij is a sum of external momenta without those from pi,... ,pj and i,... ,j is 
the default order. As the pole momentum includes at least three particles, X 23 must at least include two 
particles, X 24 must at least include one and A 25 can be empty. Similarly, X 2 ^n -3 rnust at most include 
two particles, X 2 ,n -2 must at most include one and X 2 ,n-i must be empty. Analogous restriction holds 
for all X 2 i’s in between. In this list, p 2 is the pivot momentum which is always included, and ps,... ,Pn-i 
becomes the absent momentum one by one. 

Now for (2|3] in series I, the affected two-particle poles are 


(2») 0 ex (23) 
[•3] 0 A ex [23] 


(B.3) 


where ‘ex’ means ‘except’, namely all (2») poles are filtered out except (23), and all [*3] poles are hltered 
out except [23]. The affected multi-particle pole in M is 


{P2 + X23)^ ^[2\X23\3) 


(B.4) 


since {p 2 + X 23 )'^ includes particle 2 but not 3, it is (2j3] detectable, while all other multi-particle poles are 
not, hence they remain unchanged. 

Then, for (3j4] we have 


(3*) 0 ex (34) (23) ^ (24) 

[•4] 0 A [•3] G A ex [34] 


(B.5) 


where (23) turns into (24) in H and all [*3] poles except [34] are revived, which gives a net effect that all 
[•4] poles are hltered out. Now the second multi-particle pole affected is 


(^23+^24)"^ [31P2+X2414) 


(B. 6 ) 


again, since (P 23 + -^ 24 )^ includes particle 3 but not 4, it is (3j4] detectable, while remaining multi-particle 
poles are unchanged. But don’t forget the spurious pole from (2j3], namely 


[ 2 IA 23 I 3 ) ^ [ 2 IA 24 I 4 ) or [2X24 ](x244) (B.7) 

here when A 23 includes two particles and one is 4 which is hltered out by j4), we are left with one particle 
X 24 , where Xij denotes one external momentum except that from pi,... ,pj. A nice fact is that the split 
poles [ 2 x 24 ] and (x244) are not the ones already hltered out by (2j3] and (3j4], in other words, the next 
steps will take care of them so there is no need to look back. 


- 30 - 



(B.8) 


Continuing in this fashion, for (4|5] we have 

(4.) 0 F ex (45) (34) ^ (35) (24) ^ (25) 

[•5] 0 A (no ex) 

and 

[2|X24|4) —^ [2|X25|5) or [2x2b]{x25^) 

[3|p2 + ^2414) ^ [3|p 2 + ^2515) or [32](25) (B.9) 

(^24 + ^25)'^ [4|P23 + ^25|5) 

again the split two-particle poles will be taken care of by the next steps. Note the descendent poles from 
(P 24 + -^ 25 )^, • • •, {P 2 ,n -3 + N' 2 ,n- 2 )^, P 2 n -2 longer produce split poles. After step (n — 2|n — 1], all 

descendent poles from M except split poles are 

[2|X2,n-l|n — 1) [3|p2 + -^2,n-l|n — 1) [4|P23 + A2,n-1 [^ — 1) 

(B.IO) 

[n - 4|P2,n-5 + X 2 ^n-l\n -1) [n - 3|P2,n-4 + X 2 ,n-l\n - 1) [u - 2|P2,n-3h “ 1) 
and the last step (n — Ijn] will turn them into (split poles are neglected) 

[3|Pi2|n) [4|Pi3|n) or [4|P23|n) ... 

[n - 4|Pi,„_5|n) or [n - 4|P2,n-5|n) [n - 3|Pi,n-4|n) or [n - 3|P2,n-4|n.) [n - 2|P2,n-3h) 

(B.ll) 

For two-particle poles, after series I we have 

// = {(.n)}, A= [l,...,n-l], (B.12) 

so all poles in H include |n) and any pole in A does not include jnj. Of course, all these poles can have 
orders larger than one. 

After series II, it is easy to check that 


H = {{n-l,n)], A = [1,... ,n - 2]. 


(B.13) 


For spurious poles of the form [•[ • [•) from M, series II also nicely turns them into the poles in H and A. 
To verify this, let’s single out the following series in II (note that (Ijn — 1] and (n — 2|n — 1] are set aside 
temporarily): 

(2|n-l] ... (n-3|n-l] (B.14) 

and it will be enough to take care of spurious poles [•[ • [•). Explicitly, we have 

(2|n-l]: [3|Pi2|n) [32](n - 1,n) 

(3|n-l]: [4|Pi3|n) or [4|P23|n) [43] (n - 1, n) 


(B.15) 


(n — 4|n — 1] : [n — 3|Pi^n-4|’^) or [n — 3|P2,n-4|?T-) — )• [n — 3, n — 4](n — 1, n) 
(n — 3|n — 1] : [n — 2|P2,n-3|?^) —^ [n — 2, n — 3](n — 1, n) 
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therefore these spurious poles hnally become parts of H and A after series I and II, we are only left with 
{{n — l,n)} and [1,... ,n — 2]. To concentrate them completely, series III and IV are needed, which are 
copied from I and II by replacing n by (n — 2) and swapping j*) and !•] for all deformations. 

After series III and IV, a trivial imitation gives 

H = {{n — 1, n)}, A = {[n — 3,n — 2]}. (B.16) 


Note that H is in fact completely inert to series III and IV as only A is manipulated. Therefore we manage 
to turn all poles, regardless of whether they are of two or multiple particles, physical or spurious, into a 
common denominator 

where ii,i 2 V 3)'*4 are four different arbitrary particles after a trivial relabeling. 

For n = 4, merely series I and II can turn all poles into such a form (in fact series III and IV do not 
exist). Consider the denominator there are many equivalent choices related by momentum 

conservation, after 

(2|3] (3|4] (1|3] (2|3] (B.18) 


it becomes (34)^[12]^. Again a relabeling gives the general form (iii 2 )^[^ 3 ^ 4 ]^- 

In the sequence above, series I and III are skew descendant while series II and IV are straight descen- 
dent. In general, not all physical poles can be detected by merely independent and straight descendent 
deformations. To calculate an n-particle amplitude with only these two types, one can assign a labels for 
(i| and (n — a) labels for \j] in a sequence Then two-particle poles [^ 1 ^ 2 ] with iiA 2 £ la and {jij 2 ) 

with ji,j 2 G In-a cannot be detected, where la and /„_« denote the sets of (z| and \j] respectively. 


C. (Pseudo) Polynomials and Saturated Fractions 

This part presents the classification of all three types of objects that cannot be determined by multi-step 
BCFW recursion relations. They include (pseudo) polynomials and saturated fractions. 

C.l Polynomials and pseudo polynomials 

Now, we list all (pseudo) polynomials that satisfy certain helicities, up to D = 2. This list can be similarly 
extended for D > 3. 

First consider n > 5, as there is a tricky issue of n = 4. For dimension (or polynomial degree) D = 0, 
there is only one choice of polynomial, namely 1. The helicity configuration is simply 

(0,0,0,0,0,...) (C.l) 


where ... denotes more scalars besides the minimal five. 
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For dimension D = 1, there are two choices: (••) or [••] and the helicity configuration is 





(C.2) 


For dimension D = 2, there are three choices: (••)(••), (••)[••] or [••][••]. One also needs to separate 
the cases for which one or two pair(s) of particle labels in the spinorial products are identical. Hence the 
helicity configuration can be 



1 1 

ih —, zb — , 0, . . . 
2 2 


1 1 

+ 2 >+ 2 ’ 0 ’-- 


(44±i.o,o,...) 


(± 1 ,± 1 , 0 , 0 , 0 ,...) 
( 0 , 0 , 0 , 0 , 0 ,...) 


(C.3) 


When n = 4, there exist fraction-like pseudo polynomials. Explicitly, from the following identity 


[34] 

( 12 ) 


M 

(14) 


[42] 

(13) 


M 

(42) 


[ 12 ] 

(34) 


[14] 

(32)’ 


(C.4) 


one can check that, there is no effective BCFW deformation that can detect poles above, as the numerator 
vanishes simultaneously when the denominator reaches to zero, due to momentum conservation. Moreover, 
pseudo polynomials are in fact inert to BCFW deformations, for example under (ll3]. 


[34] [34] ± z[U] 


[14] m/[U] + z 

X 


[14] [34] 


(12) (12)-z(32) (32) -(12)/(32)±z (32) (12) 

again due to momentum conservation. In general, a dimensionless pseudo polynomial takes the form 

M 

( 12 ) 


(C.5) 


(C.6) 


where x = ±1, ±2, ±3, ±4, as the spin is restricted within 2. 

Since the dimensionless pseudo polynomial behaves like a numerical factor, polynomials of n = 4 that 
share the same forms as those of n > 5 can be multiplied by these factors to generate all possible pseudo 
polynomials. Hence, for n = 4 we have the following helicity configurations: 


D = 0 : 
D = l: 
D = 2: 


{y,y,y,y) 

- 2 < y < 2 

(i 1 \ 

^ 3 

\ij^+y,2+y,y,y\ 

-2<y<- 

-y - 2 

{y,y,y,y) 

-2<y <2 

{l + y,l + y,y,y) 

-2<y<1 

/II \ 

3 ^ ^3 

[-^ + y,^ + y,y,y 

-< y < 

2 ~ 2 


(C.7) 
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where y = 0,x/2. Since (pseudo) polynomials of n = 4 are dimensionally equal, up to a global shift of all 
four helicities, for simplicity one can first select a set of representatives to analyze. As one example, the 
representative helicity configuration can be 


D = 0: 

(0,0,0,0) 


D = l: 



D = 2: 

(0,0,0,0) 

(C.8) 


(± 1 ,± 1 , 0 , 0 ) 



then y must be shifted around to cover all possible cases in (C.7). 


C.2 Saturated fractions 


The SF is an irreducible fraction, but impossible to be killed by BCFW constant extractions. A minimally 
z-power inducing constant extraction will only change its particle labels, but not its form. For example, 


[34] [56] 

( 12 ) 


(C.9) 


is a pure SF, and constant extraction (1|3] turns it into 


[34] [56] _ [14] [56] 

(12) (32) 


(C.IO) 


For a pure SF, all effective BCFW deformations are minimally ^;-power inducing. Next, consider 


[34]^[56] 

( 12 ) 


[34] X 


[34] [56] 

( 12 ) 


(C.ll) 


which is a mixed SF, namely a pure SF times a polynomial. When we use (1|3] again, it is not minimally 
^;-power inducing, i.e., 


[34]2[56] ([34]±z[14])2[56] 


0{z), 


(C.12) 


(12) (12)-z(32) 

for this mixed SF, one minimally 2 -power inducing deformation is (1|5], for which the constant extraction 
will give only one term with relabeling. But (1|3] gives more than one term, as 


[34]2[56] 

^( 113 ] 


[14] [56] [14]2(12)[56] 

(32) (32)2 


(C.13) 


note that the first term above is also a mixed SF again with relabeling. In general, under a non-minimally 
2 -power inducing constant extraction, a mixed SF will transform into more than one term. These terms 
include a special one, which is related to the original SF by relabeling. 
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To investigate SF’s more systematically, we analyze a special category: When 0 < hi,h 2 < m/2 and 
0 < /is, /i 4 , the dependence on m,Tn dissolves in (3.19), then 


D' — {—hi + Pi) + {hi + ai) + /ii + Pi) + ’y^P hj + ai), 

2 = 1,2 2 = 3,4 h<0 h>0 

and the corresponding final boundary term schematically reads 


1 




— 2hi-\-m 0i 


(12)™ [34]™ 


pi‘ n 

2 = 3,4 


h<0 h>0 


(C.14) 


(C.15) 


1st type: Let’s consider the first type of amplitudes'^, of helicity configuration {h,h,h,h,h,...) where 
0 < h < 2. It is easy to see that SF’s of n = 4 are identical to pseudo polynomials, so nontrivial SF’s only 
exist for n > 5. At this point, (C.14) becomes 


D' — {n — 4)/i -\- Pi P 2 P 'y ^ cij, (C.16) 

2=3 

and the final boundary term schematically reads 

- 47 ur(l|^M 2 |^"( 3 r 3 ( 4 |“^..(nr" x [Ij^i [ 2|^2 [ 3 | 2 / 2+«3 r 4 | 2 / 2+«4 _ _, [^| 2 / 2 +«„^ 

( 12 )^" 

where we have maximally reduced ( 1 |, ( 2 | in the fraction, hence in the numerator 

n n 

(num. of [•)) = Pi-\- P 2 + oii, (num. of [•]) = (n — 2)2h + Pi + P 2 + (C. 18 ) 

2=3 2=3 

For this type, FD arises when n = odd and h = 1/2,3/2 with D > PM arises when D' — nh = odd 

with D > and 2nh = even. From (C.16) and (C.17), it’s easy to check that: 

When D = Q, there is no consistent choice. We do not consider (pseudo) polynomials here. 

When L) = 1, n = 5, /i = 1 is SF, /i = 1/2 is FD, h = 0 is PM. 

When D = 1, n = 6 , /i = 1/2 is SF, /i = 0 is PM. 

When D = 2, n = 5, h = 1 is PM, /i = 1/2 is FD. 

When D = 2, n = 6, h = 1 is SF, /i = 1/2 is PM. 

When D = 2, n = 7, h = 1/2 is FD. 

When D = 2, n = 8, /i = 1/2 is SF. 

For D = 0,1, 2 all other cases are inconsistent since D < SF’s above are (the form of an SF is not 

necessarily unique, so we only present one example for each case) 


SF(D = l,n = 5,/i = l) 


SF(I1 = 2, n = 6, /i = 1) 


[34] [35] [45] 

( 12)2 ’ 

/MMV 

V (12) J 


SF ( D = l,n = 6,/i = - 


SF [ D = 2,n = 8,h= - 


[34] [56] 

( 12 ) ’ 

[34] [56] [78] 

( 12 ) 


(C.19) 


12 


All amplitudes in this appendix are not necessarily physical, as the investigation is purely of mathematical interest. 
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2nd type: Similarly, the second type of amplitudes is {h, h^h + 1/2, h,h ,...) where 0 < h < 3/2, then 


D' — {n — A)h + - + /3i + /32 + on. 

i=S 


For this type, we find that: 

When D = 0, there is no consistent choice. 

When L) = 1, n = 5, h = 1/2 is SF, /i = 0 is FD. 

When L) = 1, n = 6, 7, 8,..., /i = 0 is FD. 

When D = 2, n = 5, h = 3/2 is SF, /i = 1 is FD, h = 1/2 is PM, /i = 0 is FD. 

When D = 2, n = 6, h = 3/2,1,1/2, 0 is FD. 

When D = 2, n = 7, /i = 1/2 is SF, /i = 0 is FD. 

When D = 2, n = 8, 9,10,..., /i = 0 is FD. 

SF’s above are 


SF ( D = l,n = 5,/i = - j = 

SF(D = 2,n = 5,/r = ^) = [^^]W[45] 


[34] [35] 

( 12 ) ^ 

C 

(12)3 


, SF(D = 2,n = 7,/r = 2 ) = 


[34] [35] [67] 

( 12 ) 


3rd type: The third type is {h, h,h + 1/2, h + 1/2, h,...) where 0 < /i < 3/2, then 


D' — (n. — 4)/i + 1 + /?! + /32 + 'y ^ CKj. 

i=3 


For this type, we find that: 

When D = 0,1, there is no consistent choice. 

When D = 2, n = 5, h = 1 is SF, /i = 1/2 is FD, h = 0 is PM. 

When D = 2, n = 6, /i = 1/2 is SF, /i = 0 is PM. 

When D = 2, n = 7, 8, 9,..., /i = 0 is PM. 

SF’s above are 


SF(D = 2,n = 5,/i= 1) 


[34]^[35][45] 

( 12)2 


SF 


= 2, n = 6, /i 



4th type: The fourth type is {h, h,h + l,h,h ,...) where 0 < h < 1, then 

n 

D = (n — 4)/i + 1 + /?! + /32 + y ^ cni- 

i=3 


For this type, we find that: 

When D = 0, there is no consistent choice. 

When D = 1, n = 5, 6, 7,..., /i = 0 is SE. 

When D = 2, n = 5, h = 1 is SF, /i = 1/2 is FD, h = 0 is PM. 


[34] 2 [56] 
( 12 ) 


(C.20) 


(C.21) 


(C.22) 


(C.23) 


(C.24) 
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When = 2, n = 6, /i = 1/2 is SF, /i = 0 is PM. 

When = 2, n = 7, 8, 9,..., /i = 0 is PM. 

SF’s above are 

SF p = 2, n = 5, h = 1) = SF = 2, n = 6, /i = 0 


[34] [35] [36] 
( 12 ) 


(C.25) 


We will not continue to explore a fifth type, as these toy examples have provided enough intuition for the 
generation of SF’s. Naturally SF’s are highly un-physical, but their existence in some artificial (effective) 
theories cannot be excluded. Sometimes even if an SF is admitted, we can use other arguments to exclude 
it, such as its denominator contains a spurious pole. Consider an SF whose denominator is (12)^, or (12) 
but the amplitude has vanishing factorization limit under (12) —)• 0, then one needs to delicately remove 
all dependence on the spurious pole (12)^, or (12), in the known terms of (2.4) if it appears. 
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